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A search for muon neutrinos from dark matter annihilations in the Galactic Center region has 
been performed with the 40-string configuration of the IceCube Neutrino Observatory using data 
collected in 367 days of live-time starting in April 2008. The observed fluxes were consistent with 
the atmospheric background expectations. Upper limits on the self-annihilation cross- section are 
obtained for dark matter particle masses ranging from 100 GeV to lOTeV. In the case of decaying 
dark matter, lower limits on the lifetime have been determined for masses between 200 GeV and 
20 TeV. 

PACS numbers: 95.35. +d, 98.70.Sa, 96.50.S-, 96.50.Vg 



A. Introduction 

There is well-established observational evidence for the 
existence of a non-baryonic cold dark matter component 
in the universe. Although the exact nature of dark mat- 
ter remains unknown, many theories beyond the Stan- 
dard Model of particle physics accommodate stable or ex- 
tremely long-lived particles as potential dark matter can- 
didates [1 . A favored hypothesis is WIMPs (Weakly In- 
teracting Massive Particles), which are predicted to have 
masses from (9(10) GeV up to several TeV, with a theo- 
retical upper limit of 340 TeV based on unitarity require- 
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ments [2 . Particles that have properties of a WIMP nat- 
urally arise from many theories, which were not initially 
designed for this purpose. For example, the neutralino in 
supersymmetric models or the lightest Kaluza-Klein par- 
ticle in theories with extra dimensions are natural WIMP 
candidates. 

WIMPs could self-annihilate if they are Majorana par- 
ticles, or decay, resulting in stable particles like gamma- 
rays, electrons and neutrinos. A measurable fiux of such 
particles would be strong indirect evidence for the ex- 
istence of dark matter. Regions with suspected high 
density of dark matter are prominent targets for these 
indirect searches since a large dark matter density cor- 
responds to a higher flux of detectable particles. Since 
galaxies are expected to be embedded in a halo of dark 
matter, the Galactic halo and Galactic Center are exam- 
ples of such expected high density regions. 

Self-annihilation of dark matter particles in the Galaxy 
is of particular interest since a WIMP search via the 
product particles is sensitive to the velocity-averaged self- 
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annihilation cross-section {crp^v) [3]. In this sense, such 
searches are complementary to direct searches or indi- 
rect searches from the Sun. Like the direct searches, the 
latter probes the WIMP-nucleon scattering cross-section, 
since dark matter would be captured gravitationally due 
to repetitive elastic scattering. While the analysis pre- 
sented in this paper is based on a search for an excess 
neutrino flux from the direction of the Galactic Center, 
IceCube also published results from a search for dark 
matter in the Galactic halo [4 , as well as the Sun [5 . 

A search for neutrino signals from dark matter annihi- 
lations has inherent advantages over searches using pho- 
tons or charged cosmic rays since neutrinos may escape 
the production point and arrive at the detector unim- 
peded by absorption by matter or deflection by the mag- 
netic field in the central region of the Galaxy. 

B. Signal from Annihilating or Decaying Dark 
Matter 
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A description of the dark matter density in the Milky 
Way is essential for a calculation of the signal flux. There 
are various models that attempt to describe the distribu- 
tion profile of dark matter in Milky Way-size galaxies. 
These models describe a dark matter halo with a spheri- 
cally symmetric matter density, Pdm(^), that peaks at the 
center of the galaxy and decreases with distance. Profiles 
based on N-body simulations [6 of cold dark matter tend 
to show a divergent density in the central region. On the 
other hand, profiles that rely on observations of low sur- 
face brightness galaxies dominated by dark matter imply 
a flat distribution in the central region [7, 8 . Despite 
large differences in the description of the central region, 
the models show a similar density behavior at large dis- 
tances from the Galactic Center (> 8.5 kpc). 

Frequently considered models are the Navarro- Frenk- 
White (NFW) model [9 , the Moore et al. model [10^ 
and the Kravtsov et al. model [TT]. We use the NFW 
profile as a benchmark for the analysis presented in this 
paper, with the very cuspy Moore profile and the Krav- 
tosv profile, which is close to isothermal or cored profiles, 
as extreme cases. These models are compared in Fig. [l] 

The NFW profile is given by 
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with Tg = 20 kpc [9 . The normalization parameter, po? is 
chosen so that the dark matter density at the orbital dis- 
tance of the solar system {Rsc = 8.5 kpc) is the assumed 
local density Pdm(^sc) = 0.3 GeVcm"^ [3 . The choice 
of parameters is consistent with a previous analysis [4 of 
signals from the outer halo based on data from IceCube 
in the 22-string configuration (IceCube-22). 

The expected neutrino flux from self-annihilation of a 
WIMP of mass is proportional to the square of the 
dark matter density, p^M^ integrated along the line of 



FIG. 1: Dark matter density as a function of the distance 
to the Galactic Center for three halo profiles [3]. From top 
to bottom, the Moore (blue-dotted), NFW (red-dashed) and 
Kravtsov (green-solid) profiles are shown. The vertical line 
marks the orbital radius of the solar system. 
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Here ^ is the angular distance of the line of sight with re- 
spect to the Galactic Center, and Rqc and psc = p{Rsc) 
are used as scaling parameters to make Ja(^) dimension- 
less. The upper limit of the integral is 
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where i^Mw is the radius of the Milky Way dark matter 
halo. Typically, contributions from distances larger than 
the scale radius (r^) of the halo model are sufficiently 
small and can be neglected [3 . A value of i^Mw = 40 kpc 
is adopted in this analysis. 

The quantity Jao is the average value of Ja(^) for a 
chosen search region. For a circular search region defined 
by the solid angle A^l it is described by 



J. 
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Figure [2] shows Ja(^) and Jaj7 for the different halo mod- 
els and a circular search region. 

In the above scenario the neutrino flux at Earth has 
two possible components; annihilation or decay of dark 
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FIG. 2: The line of sight integral, Ja(^), and the average 
value of Ja for the solid angle defined by the opening angle, ^, 
as a function of ^ [3 . The upper x-axis indicates the fraction 
of the solid angle defined by ^ with respect to the full sphere. 
The thin lines show the value of Ja at an opening angle of ^ 
while the thick lines depict Jaq (which is the average value 
of Ja) in the cone defined by ^. The green (solid) lines show 
the Kravtsov profile, the red (dashed) lines show the NFW 
and the blue (dotted) lines corresponds to the Moore profile. 



matter. The annihilation component may be written 
as |3] 



dE 



47rm^ 



dE 
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The 1 / (47r) factor accounts for the isotropic radiation of 
the annihilation products and dA^/dE gives the differen- 
tial neutrino yield per annihilation. The factors Psc/'^x 
and Rsc normalize the dimensionless factor Ja to the 
number density, and the factor of ^ accounts for the fact 
that two particles are required for annihilation. 
The decay component is given by [12 



Rscpsc dAT^ 
47rm^ dE 



(6) 



Here r is the lifetime of the dark matter particle and 
Jd(^) is the line of sight integral over the dark matter 
density, 
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DARKSUSY \TT is used to obtain the neutrino energy 
spectra from WIMP-annihilations for various benchmark 



annihilation channels and WIMP masses. Considered are 
specific masses of 100 GeV, 200 GeV, 300 GeV, 500 GeV, 
700 GeV, lTeV,_2TeV, 5TeV and lOTeV, and 100% 
branching into 66, or r+r~ pairs. The 

obtained energy spectra represent a model-independent 
conversion of the WIMP mass to Standard Model par- 
ticles and their subsequent decay. Therefore, the same 
spectra can be used for the calculation of lifetimes in the 
case of decaying dark matter. 

Annihilation and decay to produce line spectra at 
the WIMP mass or half of the WIMP mass, respectively. 
The line spectrum is also interesting from a theoreti- 
cal point of view as it can be used to obtain model- 
independent and conservative upper bounds for dark 
matter annihilation to Standard Model final states [M]. 

The choice of the annihilation channels cover a wide 
range of neutrino spectra. The bb channel gives a softer 
energy spectrum while and u v channels provide a 

harder energy spectrum. For a specific supersymmetric 
scenario a superposition of these annihilation channels is 
expected, producing an effective neutrino energy spec- 
trum that lies between the benchmark cases considered 
in this analysis. 

The analysis presented here is based on a search for a 
muon-neutrino signal. For the calculation of the muon- 
neutrino flux at Earth all neutrino flavors are considered 
assuming a full mixing at Earth (z^g • ^^/x • z^r=l*l-l) due 
to long-baseline oscillations: 



(8) 



Here, is the initial flux at source for each flavor. 
Specific deviations from this benchmark assumption for 
vacuum oscillations, e.g. as described in [15 (^23 ^ 
7r/4,^i3 ^ 0, sin^ (26^12) = 0.86), lead to small changes 
in the predicted fluxes at the few percent level. 

Figure [3] shows the expected energy distribution of 
muon-neutrinos after full mixing at Earth for different 
annihilation channels of a 300 GeV WIMP. Annihilation 
to vv produces a line spectrum at the WIMP mass. 

The expected number of neutrino events is found by 
integrating equations ([5| or (|6| over the detector live- 
time and the direction- and energy-dependent effective 
area for annihilation or decay, respectively. This proce- 
dure permits the direct conversion of an observed flux 
at the detector to the velocity-averaged self-annihilation 
cross-section {cf^v)^ or the WIMP lifetime r. 



C. The IceCube-40 Detector 

The IceCube in-ice neutrino detector is a cubic- 
kilometer array of digital optical modules (DOMs) de- 
ployed between 1450 m and 2450 m below the surface in 
the glacial ice of the South Pole. Each DOM consists 
of a glass pressure housing enclosing a 25,4 cm diame- 
ter Hamamatsu photomultiplier tube with digitizing elec- 
tronics, and a set of LEDs for calibration p21[IZ]- The 
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FIG. 3: Energy spectra at Earth for different annihila- 
tion channels, including neutrino oscillations, for a 300 GeV 
WIMP. The spectra were obtained with DARKSUSY [l3]. 



DOMs are autonomous data collection units with a time 
resolution of less than 3 ns [18 . They detect Cherenkov 
light emitted by the medium due to interaction with rela- 
tivist ic charged particles, including charged leptons pro- 
duced in neutrino interactions in the ice, and muons pro- 
duced by cosmic-ray interactions in the atmosphere. 

The complete in-ice detector consists of 5160 DOMs 
deployed on 86 strings (60 DOMs per string). Here 78 of 
the strings are arranged in a triangular pattern and form 
a nearly hexagonal grid with 125 m inter- string spacing 
and 17 m vertical DOM separation. The remaining 8 
strings are located near the center of the detector in- 
terspersed with the standard IceCube strings, and have 
an average inter-string separation of 72 m and vertical 
DOM spacing of 7 m. Together with the surrounding 6 
IceCube strings and the central string, these strings form 
the low-energy sub-array called DeepCore [19], which is 
of particular interest for dark matter searches due to the 
lower energy threshold. 

The construction of IceCube was completed in Decem- 
ber 2010, but the detector was operational during the 
construction period from 2005 to 2010. In this paper 
we use data from the 40-string configuration of IceCube 
(IceCube-40, see Fig.[4| collected between April 2008 and 
May 2009, when IceCube did not contain any DeepCore 
strings. 

The IceCube-40 data acquisition system applied sev- 
eral trigger algorithms designed to collect different event 
topologies. This analysis considers events that satisfy the 
SMT8 trigger, a simple multiplicity trigger, requiring at 
least 8 DOMs to record a signal above threshold within 
a time window of 5 jas ^ITJ . 




FIG. 4: A sketch of string locations for the IceCube-40 de- 
tector (top view). The outer strings (open circles), which are 
used as a veto, as well as defined directions labeled as thick 
and thin are marked. 



D. Dataset and Analysis 

The data used represents 367 days of detector live- 
time [20]. Given IceCube's location at the South Pole, 
any search for a neutrino flux from the direction of the 
Galactic Center, located in the southern sky (RA 17^ 
45^ 40.04^ Dec -29° 00' 28.1''), is impacted by a sig- 
nificant background of cosmic-ray induced atmospheric 
muons and atmospheric neutrinos. The analysis strategy 
is to identify and select neutrino events which have an 
interaction vertex inside the detector volume, since at- 
mospheric muons, the dominant background, must enter 
from outside. The large instrumented volume of IceCube- 
40 allows the definition of a veto region, while retaining a 
sufficiently large fiducial volume for starting events. The 
upper thirty DOMs on each string define a top-veto while 
the complete strings in the outer layer of IceCube-40 form 
a side-veto (Fig. [4|. Events with one or more hit in the 
top veto or with the earliest hit recorded on any of the 
outer strings, are rejected. 

Events surviving the veto condition are reconstructed 
with an iterative log-likelihood track reconstruction 
method [21 . Based on the initial direction reconstruc- 
tion, only events with a zenith angle larger than 50° and 
an azimuth direction between 135° to 225° and 300° to 
30° , angular regions corresponding to the axis labeled as 
"thick" in Fig. [4j are selected. For events outside these 
regions the detector is not thick enough to define an effi- 
cient veto. Further cuts are applied to improve the angu- 
lar resolution of the data. Examples of such cuts include 
the reduced log-likelihood value from a track reconstruc- 
tion or the number of strings that have registered a hit 
in the recorded event. 

To identify events with a neutrino interaction vertex 
within the detector a likelihood algorithm has been de- 
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FIG. 5: Effective area for neutrinos from MonteCarlo simula- 
tion for the direction of the Galactic Center after application 
of the veto conditions and angular cuts (blue-dashed) and af- 
ter subsequent application of quality and signal selection cuts 
(black-solid) . 
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FIG. 6: Definition of the on-source and off-source regions. 
The extent of the declination band is 2A6 for all regions. The 
width of the on-source region in right ascension is 2Aa. The 
region between Aa and Aao is used as buffer. The region 
beyond Aao is defined as off- source. In this analysis values of 
= 8°, Aa = 9°, and Aao = 20° are used. 



veloped that analyzes the event topology [20] . This algo- 
rithm determines the likelihood of a muon to produce or 
not produce hits in the DOMs along the reconstructed 
track up-stream of an assumed starting point and it 
includes the optical properties of the ice. Maximizing 
the likelihood gives the location of the interaction ver- 
tex. The likelihood ratio of a starting track relative to a 
through-going track hypothesis is used as a classifier to 
further reduce the atmospheric muon background. 

The above-mentioned quality and signal selection cuts 
reduce the dataset by about two orders of magnitude. 
The event rate is reduced from an initial trigger rate of 
about IkHz (SMT8) to 2 Hz (initial veto and azimuth 
cuts) to 0.56 Hz (quality and signal selection cuts). Fig- 
ure [H shows the effective area of the detector before and 
after the event quality selection. The median angular 
resolution at this level is about 2.6°. The data remains 
dominated by cosmic-ray muons, undetected in the veto 
region, that enter the detector from outside. 

This analysis uses an on-source/off-source method, as 
shown in Fig. [6} to estimate the background directly from 
the data. The highest flux of WIMP annihilation prod- 



ucts is expected to arrive from the direction of a wide 
region centered at the Galactic Center. The atmospheric 
muon and neutrino background at the location of Ice- 
Cube varies with the zenith (or declination) angle but 
is approximately constant in right ascension. Therefore, 
events from a declination band around the Galactic Cen- 
ter were selected. The width of the band in declination, 
±A(5, was optimized by maximizing S/\/B^ where S and 
B are the number of expected signal and background 
events respectively for each WIMP mass and each anni- 
hilation channel. Optimal bin sizes for A^ were found to 
be between 6° and 8°. The variation of S/y/B with A5 
is rather small and a value of A5 = 8° was adopted for 
all WIMP masses and annihilation channels. 

The on-source region width in right ascension, Aa, is 
chosen to be equal to the size of the on-source region in 
declination at the zenith angle of the Galactic Center: 



cos Aa 



cos (AS) — cos^ (Ogc) 
sin^ (^Gc) ' 



(9) 



where ^gc is the zenith angle of the Galactic Center in 
local coordinates (61°). This yields a value for Aa ~ 9°. 
No sharp transition between on- and off-source regions 
exists because the dark matter density profile decreases 
continuously with distance from the Galactic Center. To 
reduce the signal contamination of the background es- 
timate from the off-source region, we define a buffer 
between the on-source region and the off-source region 
(Fig.|6|. An event is considered off-source, and thus used 
for the background estimate, if the angular distance to 
the Galactic Center is Aao ^ 20°. Events between Aao 
and Aa are discarded. 

The expected number of background events in the on- 
source region, A^bq, is: 



Aa 



180° - Aao 



(10) 



where N°^^ is the number of background events in the 
off-source region. This method of background estimation 
relies on measured data and is therefore independent of 
Monte Carlo simulations and related systematic uncer- 
tainties (e.g. modeling of optical ice properties or the 
detector response). 

Despite the daily rotation of the detector variations 
of up to 1.5% are found in the right ascension distri- 
bution of the final event selection (Fig. [t]). The origin 
of these variations is the non-uniform azimuthal accep- 
tance, caused by the asymmetric geometrical layout of 
IceCube-40, combined with periods of detector inactiv- 
ity which are not evenly distributed in time. The arising 
uncertainty is larger than the expected statistical uncer- 
tainty on the background estimate. A procedure has been 
applied to account for the uneven exposure: 

• First, we determine the normalized distribution, 
D{(j))^ of all measured event azimuth angles, 0, 
in detector coordinates that do not originate in 
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FIG. 7: Relative variation of the event distribution in RA 
(with respect to the Galactic Center, thus the x-axis is la- 
beled ARA). N is the number of measured events per bin, 
(N) is the average number of events in all bins. The mea- 
sured (black-dotted) and expected (red-dashed-dotted) distri- 
butions are shown, along with the deviation from expectation 
(blue). 



TABLE I: Jaq for annihilating dark matter for the considered 
halo density profiles and the chosen on-source region. 





NFW 


Kravtsov 


Moore 


Jaq 


240 


23 


900 



spends to a rotation of the azimuth pdf by an angle, 
A(j)i^ that can be calculated from the event time. 

• Lastly, we calculate the sum of all the obtained 
distributions (in RA): 



N 



N 



D'{RA) = Y,D[ = Y,D{RA + ^4>i), (H) 

i i 

where N is the total number of events. 

The application of this procedure results in an exposure- 
corrected expected distribution of background events in 
right ascension, normalized to the number of events in the 
full dataset. Figure [7| shows the accuracy of the number 
of predicted events. The agreement between the number 
of measured and expected events is at the 0.1%-level. 
This approach also takes seasonal variations into account. 
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FIG. 8: Initial azimuth pdf containing all events from the off- 
source region. The cuts on the reconstructed azimuth were 
performed on an initial reconstruction. This figure shows the 
azimuth distribution based on an advanced reconstruction. 
The cut regions excluding events from the "thin" (Fig. [4| 
detector orientation are clearly visible. 



the on-source region. This distribution, shown in 
Fig. |8j corresponds to the effective probability den- 
sity function (pdf) for the observation of an event 
from a certain detector direction, integrated over 
the exposure in time. This pdf rotates with the 
detector and is not the exposure in right ascension. 

• Second, we transform for each event time the full 
pdf to equatorial coordinates (RA,(^): D'-{RA) = 
D{RA -\- A(j)i). At IceCube's location, this corre- 



E. Systematic Uncertainties 

Predictions of the expected signal are dominated by 
large variations in the description of the central region 
of halo models (Fig. [T]). In particular, in the direction of 
the Galactic Center, the line of sight integral, J(^), spans 
several orders of magnitude (Fig. [2|. For this reason we 
have chosen the NFW profile as a benchmark model, and 
do not consider the halo-model dependence a systematic 
uncertainty. 

The limits obtained in this analysis can be rescaled to 
limits for the Moore or Kravtsov profile by multiplication 
of the limit for each mass or channel with Jaq^/'^aq^ 
where J^^^ is the desired profile. The corresponding 
values of J/\q were obtained for the actual box-shaped 
search bin (not by means of equation Q) and are shown 
in Table [Tj However, especially for flat-cored profiles the 
signal-contamination of the off-source region has to be 
considered when rescaling the limits. 

The background prediction is derived from the high- 
statistics off-source data sample. The systematic uncer- 
tainty in this prediction is limited to small directional 
variations in the cosmic ray flux, e.g. as reported in [22]. 
A Kolmogorov-Smirnov test for compatibility of the rel- 
ative deviations to a normal distribution yields a p- value 
of 0.77. Within the limited statistics of 72 values (5°- 
bins) , the result does not hint at a possible anisotropy in 
this data sample. 

Though the background estimate is data-driven, sys- 
tematic uncertainties arise from the Monte Carlo-based 
estimation of the detector sensitivity for signal events. 
Typically, the energies considered in this analysis are 
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TABLE II: Summary of systematic effects. 





Contribution 


live-time 


< 1% 


v cross-section and 




propagation 


7% 


photon propagation and 
optical ice properties 


20 %-25 % 


total (conservative) 


33% 



lower than for other IceCube analyses, e.g. the search 
for high energy extr agalactia neutrinos, and are largely 
similar to the energies considered in Ref. [U 0. For 
lower energies, the uncertainties on the neutrino cross- 
section and neutrino and muon propagation are on the 
order of 7% [5^. In addition, uncertainties corresponding 
to the detector up-time are less than 1%. Uncertainties 
related to the detection efficiency of low energy events, 
such as the sensitivity of optical modules [16 and the 
photon propagation through the ice [23 , are consider- 
ably larger. These effects lead to a typical uncertainty 
on the total effective area of about 25 % for lower en- 
ergies, (9(100 GeV), and about 20% for energies near 
1 TeV. Though the impact of the above-mentioned sys- 
tematics is energy-dependent we conservatively assume a 
33% effect on the signal detection efficiency, edet- Here, 
edet corresponds to the integration over energy of the 
neutrino multiplicity per annihilation multiplied by the 
effective area: 



/ 



dE 



(12) 



The value of 33 % is propagated into the limits by means 
of equation ([5|, 



{(JAV) (X 



lim 



edet ■ (1 ± 0.33) ' 



(13) 



for all channels and masses, where A/'nm is the limit on the 
number of signal neutrinos in the sample. The systematic 
uncertainty on the lifetimes have been included by means 
of equation ([6|, 



r (X 



edet • (1 ± 0.33) 



lim 



(14) 



The presented limits have been shifted upwards by a fac- 
tor of 1/0.67 to include systematic effects in a conserva- 
tive way. 



F. Results 

The number of observed events in the on-source re- 
gion was 7V°" = 798 842, while N^l = 798 819 back- 
ground events were expected based on estimates from 
the approximately 5 times larger larger off-source region, 
using equation (10). This procedure yields a relatively 




[GeV] 



FIG. 9: Upper limits (90% C.L.) on the self-annihilation 
cross-section for the four channels at nine WIMP masses. The 
hashed area represents the natural scale at which WIMPs 
are thermal relics [25]. The black-dotted line is the unitarity 
bound [21 [26]. The four limits were derived assuming the 
NFW density profile. 



small uncertainty on the background expectation. If the 
null hypothesis is valid, the outcome of a large number 
of measurements of N°'' is expected to follow a Pois- 
son distribution around N^^- The small difference of 
A7V = N"""" — A/'°J — 23 events, compared to a stan- 
dard deviation of about 916 events, is compatible with 
no neutrino signal from annihilations of dark matter par- 
ticles. Based on Feldman- Cousins [2T, an upper limit 
(90% C.L.) of 1504 signal neutrinos is computed. From 
equation ([5|, the limit on the velocity-averaged WIMP 
self-annihilation cross-section {(tav) for annihilations in 
the Galactic Genter is derived for various WIMP masses 
and annihilation channels (Fig. |9|. Similarly, with equa- 
tion ( |6l), lower limits on the WIMP lifetime are calculated 
(Fig.™. 



G. Comparison with Results from Other Indirect 
Searches 



Several satellite-borne and balloon experiments have 
reported deviations of measured cosmic-lepton spectra 
from theoretical predictions. The PAMELA collabo- 
ration observed an increased positron fraction at ener- 
gies above ~ 10 GeV [28], that disagrees with the con- 
ventional model for cosmic-rays [29 (secondary positron 
production only). This measurement has also been con- 
firmed by the Fermi-LAT collaboration ^30j. Further, 
the Fermi-LAT [31] and HESS [32l [33] measurements 
of primary electrons showed a discrepancy when com- 
pared to conventional diffuse models. Though the results 
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[GeV] 



cross-section and the muon range [38] . 

In a different approach, the dark matter self- 
annihilation cross-section can also be constrained based 
on precise measurement of the cosmic microwave back- 
ground (CMB). A too large self-annihilation cross-section 
can cause significant distortions of the CMB depending 
on the annihilation mode. The lack of such distortions 
can be translated to a limit on (aAv). 

Very stringent limits were obtained [39l [40] from mea- 
surements of temperature, anisotropy and polarization of 
the CMB, based on WMAP [27 and ACT [41] data. 

We refrain from comparing the derived limits to our 
results as our analysis targets dark matter properties in 
the present epoch, while conditions in the early universe 
were substantially different, e.g. the velocity distribu- 
tion, as well as the energy absorption processes in the 
photon-baryon plasma [42] . 



FIG. 10: Lower limits (90% C.L.) on the lifetime for the four 
channels at nine WIMP masses. The horizontal line repre- 
sents the age of the universe [27 . The four limits were derived 
assuming the NFW density profile. 



may be interpreted as signals from nearby astrophysical 
sources they could also be explained by annihila- 
tion or decay of leptophilic dark matter with a mass of 
0{1 TeV). With this assumption, preferred regions in the 
{aAv)-m^-pldine can be identified [35^, ^ and compared 
to the limits obtained in this analysis. Figure 11 shows 
channel limits compared to those regions. Fig- 
12] shows a similar comparison for the r+r~ channel 



X X -> 



our 



ure 

where the current limits are close to the preferred regions 
of the combined PAMELA and Fermi data. The better 
exclusion power in the r-channel is not due to a more 
stringent limit, but rather the location of the favored re- 
gions since there is significant energy loss during electron 
propagation, any electron or positron signal must orig- 
inate from nearby sources ((9(1 kpc)). Given this, the 
shape of the preferred regions does not depend strongly 
on the halo profile. 

Recently, the Fermi-LAT collaboration [37] reported 
results from a search for dark matter annihilation sig- 
nals from 10 dwarf galaxies. For gamma-ray experiments 
dwarf galaxies are very promising targets due to a rela- 
tively high signal-to-noise ratio. The obtained limits are 
very stringent and a comparison to our Galactic Center 
limits, as well as the IceCube-22 Galactic halo limits [4j, 
is shown in Fig. 
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At lower WIMP masses the dwarf 
limits are more constraining than the IceCube limits by 
several orders of magnitude. However, the IceCube lim- 
its on direct annihilation to neutrinos can be regarded as 
a very conservative bound on the total dark matter an- 
nihilation to Standard Model final states [T4^. At higher 
WIMP masses (> TeV), where there is larger neutrino 
cross-section and muon range, IceCube benefits from the 
large instrumented volume due to the increasing neutrino 
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FIG. 11: Upper limits (90% C.L.) on the self-annihilation 
cross-section for the annihilation channel to muons at nine 
WIMP masses assuming the NFW density profile (solid line) . 
The hashed area represents the natural scale at which WIMPs 
are thermal relics . The black-dotted line is the unitarity 
bound [21126]. The green (3cr) region is based on the PAMELA 
excess, the red region represents the PAMELA data combined 
with the Fermi and H.E.S.S data (3cr and 5cr) [35j . 



H. Summary and Outlook 

The analysis of IceCube data in the 40-string config- 
uration has been presented and yields limits for a wide 
range of WIMP masses and annihilation channels. Fig- 
ure 



13 shows the limits from this analysis compared to 



those from an analysis of the outer halo with IceCube- 
22 [4 . Roughly an order of magnitude improvement has 
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FIG. 12: Upper limits (90% C.L.) on the self-annihilation 
cross-section for the annihilation channel to tau leptons at 
nine WIMP masses assuming the NFW density profile (solid 
line). The hashed area represents the natural scale at which 
WIMPs are thermal relics [25]. The black-dotted line is the 
unitarity bound [21 [26] . The green (3cr) region is based on the 
PAMELA excess, the red region represents the PAMELA data 
combined with the Fermi and H.E.S.S data (3cr and 5cr) [35J. 

been achieved in the lower mass range, while at higher 
masses neutrino limits provide some of the most strin- 
gent constraints. Further, neutrino constraints can be 
used to test scenarios motivated by the anomalous lep- 
ton observations and are also extremely effective in ruling 
out scenarios for high mass WIMPs. Recent Fermi ob- 
servations of diffuse gamma-ray flux could hint at such 
scenarios, in which the unitarity bound does not apply, 
for example in p-wave dominated annihilation scenarios, 
or models of composite dark matter [43 . 

The sensitivity of this analysis is limited by the signif- 
icant atmospheric muon background in direction of the 
Galactic Center as well as the relatively high detection 
threshold. Substantial improvement is expected from 
analyses using the considerably larger detector config- 
uration with 79 strings, IceCube-79, and the complete 
IceCube-86 detector. A larger detector provides a larger 
fiducial volume while maintaining a sufficiently large and 
effective veto volume to suppress the atmospheric back- 
ground [44 . In addition, preliminary estimates suggest 
that the DeepCore low-energy extension, first present in 
IceCube-79, increases the effective area significantly in 
the energy region between 10 GeV and 100 GeV. Dedi- 
cated data filters for the Galactic Center have been devel- 
oped and implemented for IceCube-79 [38 , and further 
improved in terms of signal acceptance for the complete 
IceCube-86 detector. 



FIG. 13: The upper limits from this analysis compared 
to those from an analysis of the outer Galactic halo using 
IceCube-22 data i4j, as well as limits obtained from observa- 
tion of dwarf galaxies by the Fermi-LAT collaboration [37] . 
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